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Abstract. We introduce the (logarithmic) bipartite fidelity of a quantum system 
A U B as the (logarithm of the) overlap between its ground-state wave function and 
the ground-state one would obtain if the interactions between two complementary 
subsystems A and B were switched off. We argue that it should typically satisfy an 
area law in dimension d > 1. In the case of one-dimensional quantum critical points 
(QCP) we find that it admits a universal scaling form ~ hi£, where t is the typical 
size of the smaller subsystem. The prefactor is proportional to the central charge c 
and depends on the geometry. We also argue that this quantity can be useful to locate 
quantum phase transitions, allows for a reliable determination of the central charge, 
and in general exhibits various properties that are similar to the entanglement entropy. 
Like the entanglement entropy, it contains subleading universal terms in the case of a 
2D conformal QCP. 
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Introduction. — A major challenge in the study of quantum many-body systems in 
condensed matter physics is the understanding and characterization of new exotic phases 
of matter, such as quantum critical or topological phases. For this purpose, various 
quantities and concepts have been introduced, some coming from quantum information 
theory. Amongst them, one of the most heavily studied is the entanglement entropy 
(EE)p], defined through a bipartition of a total system A U B, usually in a pure state 

S = -Trp A \np A , p A = Tr B (1) 

The EE of a ground state is known to be universal at one-dimensional quantum critical 
point (QCP) [21 [31 H], and the leading term allows for an accurate determination of 
the central charge. In higher dimension d > 1, it obeys an area law (with possible 
logarithmic corrections |5J): if L is the typical size of the smaller subsystem, then S 
scales as L d ~ l . In this case, subleading terms[6j [91 [10] encode universal features of 
the system, characterizing quantum criticality or topological order. Despite all of these 
theoretical works, connecting the EE to experimentally measurable quantities remains a 
formidable task. One reason for that is that characterization of entanglement in a many- 
particle system generically requires the measurement of a prohibitively large number of 
observables. 

A different class of quantities is the one of overlaps, or fidelities. The idea is perhaps 
more intuitive and can be traced back to Anderson's orthogonality catastrophe [II] . Let 
H{\) be a Hamiltonian which depends on a physical parameter A that can be varied. If 
| A) is its ground state, the fidelity is 

/(A,A') = |(A|A')| 2 . (2) 

Close to a QCP the fidelity susceptibility x(A) = (dyf(X, A')) A , =A diverges, so this 
quantity can be used to detect quantum phase transitions [T2], [131 EH]- The scaling 
behavior of the fidelity has been studied in various systems both analytically and 
numerically [151 US]- Overlaps are also interesting when considering time-evolutions. 
Starting from an initial state 1^(0)), one can ask what the overlap of the wavefunction 
with the initial state after time t is. The Loschmidt echo C{t) = \ (ip(0)\ip(t)) | 2 [171 [18] 
has been studied in connection with NMR experiments [19], and also in the context of 
quantum criticality [201 ETJ • 

In this letter, we introduce an overlap which shares some common properties with 
the EE (in particular we keep the idea of cutting the system into two parts), despite 
being conceptually simpler. We call it logarithmic bipartite fidelity (LBF), and claim 
that it provides valuable insights into quantum critical phenomena. We shall see in 
particular that the LBF obeys an area law in d > 1, can be useful to locate QCPs, and 
has a universal scaling form at one- dimensional QCPs that involves the central charge 
c, much in the spirit of the EE. 

Bipartite fidelity. — Let us consider an extended quantum system AU B described 
by the Hamiltonian 

H = H A + H B + H^ub (3) 
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Figure 1. Bipartition of a 2D system (on the left), along with the three partition 
functions of Eq. [5] in d + 1=31?. Region A is in blue when decoupled from B. 



where [Ha, Hb] = and H A ^ B contains all the interaction between A and B. We denote 
by \A) (resp. \B)) the ground-state of Ha (resp. Hb), by |A®.B) = |A)<g)|-B} the ground- 
state of Ha + Hb, and by \A U B) the ground-state of H. We introduce the bipartite 
fidelity \ (Al} B\A® B)\ 2 , the overlap between the ground-state of the total Hamiltonian 
H, and the ground state of a Hamiltonian Ha + -ffs where all interactions between A 
and B have been switched off. A more physical way of looking at this quantity is to 
interpret it as a probability of measuring a given energy after a local quantum quench. 
Let us imagine that the system is initially disconnected (i.e. it is in the ground state of 
Ha + Hb), and that at time t = the interaction between A and B is instantaneously 
switched on. Then, at time t > 0, the system evolves with the total Hamiltonian ([3]). If 
one measures the energy of the system just after the quench, the probability of finding 
the ground state energy is given by | (AU B\A ® B) | 2 , which is the bipartite fidelity. For 
later convenience, we consider (minus) the logarithm of this quantity 

F A ,B = -ln(\(AUB\A®B)\ 2 ), (4) 

and call it logarithmic bipartite fidelity (LBF). The symbol J^a,b is chosen because it 
can be interpreted as a free energy in a classical system. 

Free energy and area law. — The LBF is nothing but a linear combination of free 
energies of different d + 1-dimensional systems. Indeed, in a euclidean picture, the 
ground state |0) of a Hamiltonian H can be seen as the result of an infinite (imaginary) 
time evolution starting from any state \s), provided (0|s) ^ 

e~ TH \s) ~ e- rEo \0}{0\s}. (5) 

T— > + 0O 

E is the ground state energy. Making use of this, our scalar product can be expressed 
as a ratio of classical d + 1-dimensional partition functions 

(A I I B\A ® B) = lim —— Za > b ^ ( 6 ) 
r^oo ^Z Aub (t)Z a ® b (t) 

Za®b — (s \ e~ 2T ( HA+HB ^ \s) is the partition function corresponding to two independent 
systems A and B, Zaub — (s \ e~ 2T ^ HA+HB+H A^B) \s) is the partition function of the total 
system, and Za,b = (s\ e~ T<yHA+Hl3 ^ e~ T<yHA+HB+HAi ~> B ^ \s) corresponds to the case when A 
and B are decoupled from — r to 0, and coupled afterwards (see the d = 2 example in 
Fig. [1]). In terms of free energies / = — In Z, the LBF is then 
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J 7 A,B — 2fA,B — f A®B — IauB- (7) 

The different terms in ([7]) are expected to be extensive in the thermodynamic limit. 
There is a bulk free energy fd+\ per unit volume, a "surface" free energy fa, and a 
"line" free energy fd-i- The bulk and surface energies are canceled out by the linear 
combination, and we get 

T A ,B = f d -iL d - 1 +o(L d - 1 ), (8) 

where L d ~ l is the "area" of the boundary between A and B in the initial c?-dimensional 
system. This is the area law for the LBF. We expect this to be true for generic systems, 
as is the case usually for the EE, however like for the EE [5], exceptions are possible. 

ID conformal QCPs. — In general, Fa,b should be finite away from a QCP, because 
the correlation length £ is small and the two ground-states are very close to each other, 
except on a thin region of typical size £. At a QCP however, this is no longer true and 
Ta,b can become large. As the calculation of the overlap boils down to a free energy, the 
scaling behavior should be controlled by the conformal symmetry only. This is indeed 
the case, and this result constitutes the central point of our work. The geometries 
considered are shown in Table, dj along with the corresponding formulae we derived for 
the bipartite fidelity. For example, geometry (a) consists in a finite chain of length £ 
connected to another finite chain of length L — i. When £ -C L, the scaling of the fidelity 
takes the following simple form 

J- AjB ~£ln£ (9) 

o 

This result is similar to the one for the EE [3]: 5* ~ (c/6) ln£. Another simple result is 
for the symmetric case £ — L/2, where 

F A , B ^lhxL y (10) 

o 

whereas the EE behaves as S ~ (c/6) InL [IJ. There is no general relation between the 
LBF and the EE though. Our analytical results (Table. [I]) do not match the ones for the 
EE S ~ |ln[^sin^], which can be traced back to the fact that the Cardy-Calabrese 
derivation [3J of the EE involves a local twist operator, whereas the LBF cannot be 
expressed as a correlator of a local field. Both quantities have a similar qualitative 
behavior however. 

Conformal field theory derivation. — The results in Tab.[T]assume that the boundary 
conditions at the boundaries of A U B (if any), A and B are conformal boundary 
conditions [22] • Moreover, for simplicity, we assume that these boundary conditions 
are the same everywhere. For different boundary conditions, the scaling dimensions of 
the different boundary condition changing operators [22J would modify our results [23J . 
With those two assumptions at hand, the two special cases of geometry (a) given in 
(19~|)-f lTU|) are straightforward applications of the celebrated Cardy-Peschel formula [24J: 
in the three geometries shown in Fig. [2] there is one corner with angle 2ir and several 
corners with angle at infinity. The contributions at infinity cancel, and we are left 
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InL + g b (x) + g b (l - x) 
3 - 6x + Ax 2 



6(1 -x) 
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Table 1. Geometries considered, along with the leading term for the two LBFs (T a 
and J-f,), as a function of L and x = t/L. 
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Figure 2. The three geometries giving the terms Ja,b, Jaqb and Jaub for the case 
(a) in Tab. [U The blue intervals [—A, +A] are used for the regularization. 



with the contribution of the corner with angle 2ir, which gives a total factor of 2 x ^ In I 
(or 2 x ^lnL). 

In the general case (a) in Tab. [lj the calculation goes as follows. /a,b, Ia®b and 
Jaub in (El) are the free energies in the geometries shown in Fig. |2j Let w = x + iy 
be the complex coordinate such that the lower (resp. upper) boundary corresponds 
to y = Qmw = (resp. y = L). Let us consider the mapping w h-» w + i5£ if 
Smw G (0, L), and w i-> w otherwise. Such a mapping keeps L fixed but changes i into 
t + <W. The variation of the free energy can be expressed in terms of the T yy component 
of the stress-tensor as [22] 



0Ja,b = hm — 

A— s-oo Z7T 



dx, (11) 



-A 

and there are similar expressions for /a®b and /aub- Each of these expressions diverges 
when A — > oo, but the combination ([7j) is finite. To evaluate <5/a,b one needs the stress- 
tensor in the pants-like geometry (Fig. [2] left) . We use a conformal mapping z i— > w(z) 
from the upper half-plane z6C, z > to the pants-like geometry 

£ L— I ( I \ 
w(z) = - ln(l + z) + In - 1 . (12) 

7T 7T \ L — £ J 

In the half-plane, one has (T(z)) = 0, so using the transformation law of the stress- 
tensor [22] we get (T(w)) = --^(dw/dz)~ 2 {w, z}, where {w, z) = ^ — § (^) is the 
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Schwarzian derivative of w with respect to z. Since (T yy ) = (T(w)) + (T(w)}, we find 

where w(x\) = A, w{x2) = —A, w(xz) = A + iL, w(x^) = —A + iL. In a strip of width 
i the stress-tensor is [22] (T(w)) = so 5f A ® B = 2A^ff (1/(L - £) 2 - l/^ 2 ) and 

SJaub = 0. Finally, introducing the aspect ratio x = £/L and taking the A — > oo limit 
in the sum (J7J), we get the following equation for J 7 = Ta,b '■ 



5 J 7 _ c 
5x 6 



x 2 (2 —x), x 2 — 1 , . s 1 x 
— i L In x -| __ l n (l - x) H 



— x -)■ 1 



2(1 -x) 2 2x x 7 41 -x_ 

This can be integrated to give the formula (a) in Tab. HJ For the periodic case (b), 
we need a conformal transformation which maps the upper half-plane onto a cylinder 
with two slits 

w(z) = ±\n(l + ^(z 2 -lf)-j\nz\ (14) 

and the result follows from a similar calculation. 

Numerical checks. — We study the XY-chain in transverse field, with boundary 
conditions o~ x L+1 = o~ v L+l = 0: 

H =-jl (^^^+1 + ^y^M+i + haij • (15) 

Two cases are of special interest: r = and h = is the critical XX chain, in the 
universality class of the free boson (c = 1); whereas r = 1 is the Ising chain in transverse 
field (ICTF), critical at h — 1 with c = 1/2. Using a Jordan Wigner transformation, 
Ha, Hb and Haub can be recast as free fermions Hamiltonians, and diagonalized by 
a Bogoliubov transformation. Keeping track of the changes of basis, the overlap can 
be expressed as a fermionic correlator, and reduced to a L x L determinant after some 
algebra. 

Results for the XX chain are shown in Fig. [3] for geometry (a), and agree very well 
with the CFT prediction. We also checked our formula for geometry (b). 

Results for the overlap as a function of h in the ICTF are also shown in Fig. [3j 
The quantum phase transition at h = 1 can clearly be seen, even with relatively 
small system sizes. In the vicinity of the critical point, the correlation length is 
known to diverge as £ ~ \h — 1| _1 , and the rescaled overlaps can be made to collapse 
onto a universal curve. Similar to the EE [25], we also have the exact relation 
T {XX) {LJ) = 2J r ( Isin §)(L/2,£/2, h = 1) on the lattice. 

Time evolution after a local quench. — Let us consider again the system in Tab. Q] 
(geometry a). It is prepared at time t = in the state \A <g> B). Then for t > the two 
parts A and B interact, and the system AUB evolves with the total Hamiltonian ([3]). It is 
well-known that the EE for such a system grows as [26] S ~ | In t for a Vpt £, L (a 
is the lattice spacing and vp the Fermi velocity). This logarithmic growth has given rise 
to speculations about a possible relation between the EE and the statistics of fluctuations 
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Figure 3. Left Panel: XX chain numerical results for the LBF T = Ta.b in 
geometry (a). Right Panel: Geometry (a) with i — L/2 for the ICTF. Rescaled 
LBF FyLQi — 1)) — (c/8) lnL can be seen to collapse onto a single universal curve in 
the vicinity of the critical point. Inset: J-(L, h) as a function of h. 



of the current between the two parts in certain fermionic systems |27J, which would open 
the route to an experimental measure of the EE. Here we stress the fact that the time- 
dependent LBF grows logarithmically as well. Actually, the bipartite fidelity in that 
case is nothing but a Loschmidt echo 

C(t) = \(A®B\ e iHt \A ® B) | 2 (16) 

so the LBF is J-A,B(f) — — ln£(£). Its universal behavior can be derived in CFT as 
follows. In imaginary time the scalar product (A <g> B\ e~ rH \A <g> B) is the partition 
function of a 2D statistical system in a strip with two slits separated by a distance VpT. 
In the limit vpr <C £, L the two slits almost touch each other. Again, the Cardy-Peschel 
formula shows that the contribution to the free energy of each of these corners scales 
as £: In Vft. The LBF behaves then as J z a,b{'t) ~ § hi|r|. Going back to real time 
T — y 6 — it, we find 

^AB(*)~|ln(l + 5) ~ i ln *" ( 17 ) 

2D conformal QCPs. — As discussed before, in 2D the bipartite fidelity should scale 
linearly with the system size J 7 = j\L + o(L), where f\ depends on the microscopic 
details of the theory. Universal quantities, if there, have to be looked for in subleading 
corrections. We consider the simple example of a critical quantum dimer wave functions, 
whose amplitudes are given by the Boltzmann weights of a 2D classical dimer model. 
In the continuum limit this wave function is related to a free boson CFT with 
compactification radius R [29 , 28 j • For the geometry of a cylinder of height L y ^> L x 
cut into two parts (see Ref. [10]), the LBF can be expressed using classical partition 
functions for the dimers [ID], T = In Z DD (L X , L y ) - 2 In Z DD (L X , L y /2). D stands for 
Dirichlet and encode the conformal invariant boundary condition at both end of the 
cylinders in the continuum limit. The first subleading term is a constant related to the 
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(Dirichlet) Affieck-Ludwig boundary entropy |30j sp computed in Ref. [51] : 

JF~ f x L-2s D = hL + \nR. (18) 

It would certainly be interesting to study this idea in more complicated models. 

As is the case for the EE [6], we also speculate that subleading terms in the LBF 
might be used to identify topological order. For certain trial wave-functions, such as 
Rokhsar-Kivelson triangular lattice quantum dimer and Levin- Wen string net wave- 
functions, the LBF is nothing but the n — > oo Renyi entropy. Then the arguments 
in [6] would yield the same subleading constant in the LBF and in the EE. Another 
way of looking at this would be to compare the eigenstate of the reduced density 
matrix associated with its largest eigenvalue to the actual ground state of the physical 
Hamiltonian Ha ■ For a generic Hamiltonian there is no reason why there should be any 
relation between those two states. However, for special Hamiltonians associated with 
trial wave functions, we speculate that they might be closely related to each other. We 
leave this important open question for future studies. 

Conclusion. — We have introduced the LBF of an extended quantum system AUB, 
and studied some of its properties. We have shown in particular that it generically 
obeys an area law, and exhibits universal behavior at ID and 2D QCPs, like the EE. 
We note that its simple definition makes it easier to grasp intuitively than the EE, and 
convenient to study using standard analytical and numerical methods. This could be 
particularly useful in dimension d > 1, where Quantum Monte Carlo algorithms allow 
to compute efficiently ground state overlaps, as opposed to the von Neumann entropy. 
Therefore, we believe that the LBF can be a useful and general tool in the study of 
quantum many-body systems. 
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L. Jacobsen, G. Misguich, V. Pasquier and H. Saleur for valuable discussions and 
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